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Abstract: The holographic approach to building the p-wave superconductors results in 
three different models: the Maxwell-vector, the SU(2) Yang-Mills and the helical. In the 
probe limit approximation, we analytically examine the properties of the first two models 
in the theory with dark matter sector. It turns out that the effect of dark matter on the 
Maxwell-vector p-wave model is the same as on the s-wave superconductor studied earlier. 
For the non-Abelian model we study the phase transitions between p-wave holographic 
insulator/superconductor and metal/superconductor. Studies of marginally stable modes 
in the theory under consideration allow us to determine features of p-wave holographic 
droplet in a constant magnetic field. The dependence of the superconducting transition 
temperature on the coupling constant a to the dark matter sector is affected by the dark 
matter density pn. For po > p the transition temperature is a decreasing function of a. 
The critical chemical potential pc for the quantum phase transition between insulator and 
metal depends on the chemical potential of dark matter pp) and for = 0 is a decreasing 
function of a. 

Keywords: Gauge-gravity correspondence, Holography and condensed matter physics 
(AdS/CMT), Black Holes 


^rogat@kft.umcs.lublin.pl, marek.rogatko@poczta.umcs.lublin.pl 
^ karol@ ty t an .umcs.lublin.pl 




Contents 


1 Introduction 1 

2 Models of p-wave holographic superconductor with dark matter sector 3 

2.1 Model of SU(2) Yang-Mills p-wave holographic superconductor with dark 

matter sector 5 

3 P - wave holographic metal/superconductor phase transition 6 

3.1 Condensation value 9 

4 Insulator/ holographic p-wave superconductor phase transition 11 

4.1 Critical chemical potential 12 

4.2 Critical phenomena 14 

5 P-wave holographic droplet 16 

6 Summary and conclusions 18 


1 Introduction 

The gauge/gravity duality provides a powerful theoretical method which enables a better 
understanding of the strongly coupled systems [l]-[3]. Originally proposed as the equiv¬ 
alence between type IIB superstring theory on AdS^ X spacetime and AA = 4 SU{N) 
supersymmetric Yang Mills theory on (3 -|- l)-dimensional boundary, has later been gener¬ 
alized to other gravitational backgrounds [4]. The correspondence empowers an equality 
between the quantum field theory in d -dimensional spacetime and the gravity theory in 
(d -|- l)-dimensions, and its usefulness originates in a strong-weak duality [5]. Namely, the 
gravity dual of the strongly coupled quantum field theory is tractable in a perturbative 
approach. 

The AdS/CFT correspondence has been recently proposed as a method to describe 
superconducting phase transition of the single s-wave superconductor [6] . Shortly afterwords 
it has been generalized to take into account other symmetries, like simple p or d wave. 
The generalizations require proper choice of the condensing field and appropriate gravity 
background. To describe d-wave holographic superconductor the charged massive spin-two 
held in the bulk [7]-[10] is required. The works related to building the holographic p-wave 
superconductor [11] have indicated a number of equally feasible possibilities and resulted 
in the multitude of approaches to study them. Inter alia the hve-dimensional supergravity 
framework [13], and the Sturm-Liouiville eigenvalue problem [14] have been applied. A 
handful of novel results [12, 15, 16] have also been reported. Especially intriguing result [15] 
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is the change of the order of the superconducting transition showing up when backreaction 
is taken into account. The second order phase transition is replaced by the first order one, 
when matter field couplings are beyond a critical value. 

Recently, the aforementioned studies were generalized in many other ways. The modi¬ 
fication of gravity theory by considering the five-dimensional AdS solitonic metric has been 
proposed [17]. It enabled the construction of the holographic insulator/superconductor 
phase transition at zero temperature [18]. Namely, the AdS soliton line element dual to a 
confined field theory with a mass gap, imitates an insulator phase [19]. The strength of 
various kinds of matter backreactions has been shown to generate new phase transitions 
[ 20 , 21 ]. 

The marginally stable modes of scalar and vector perturbations of the AdS spacetime 
have revealed the outset of the phase transition and help to study the influence of magnetic 
field on them [22, 23]. In agreement with known phenomenology the magnetic field in the 
holographic theory makes the phase transition harder to occur. Both analytical and nu¬ 
merical methods devoted to the properties of s-wave and p-wave insulator/superconductor 
phase transitions were investigated in [24, 25]. The studies in Gauss-Bonnet gravity were 
presented in [26, 30], the effects of the Weyl corrections on p-wave holographic phase transi¬ 
tions [27] were studied in [28]. The p-wave holographic superconductors in different gravity 
backgrounds [29], and in the presence of non-linear electrodynamics [31, 32] and other 
non-trivial conditions [33-35] have been elaborated. 

The important problem being of special interest for the present work is a possible 
matter configuration in the Universe. According to numerous studies more than 24% of the 
matter is invisible and therefore dubbed as dark matter. There exist various proposals of 
how to model this component of matter. In this paper the point of view has been accepted 
according to which the dark matter is described by the U(l) field [21] analogous to the 
Maxwell one and coupled to the ordinary matter. The coupling constant a is treated as a 
free parameter, which value is bounded by | a |< 2. 

This model of dark matter is supported by numerous astrophysical observations [36- 
40] and other experimental data related to the muon anomalous magnetic moment [41] and 
experimental searching for the ’dark photon’ [42-46]. 

The main aim of this paper is to study the influence of the dark matter sector on 
the properties of holographic p-wave phase transitions. These studies may hopefully result 
in discovery of some qualitatively untypical behavior which could be tested in the future 
experiments elucidating the cloven nature of dark matter. 

The superconducting transition is signaled by spontaneous breaking of the U (l)-gauge 
symmetry. In the case under consideration the rotational symmetry is also broken by a 
special direction of a vector field, which is obtained by the condensation of a charged vector 
field. The present paper is the generalization of the previous works [47-49], were various 
aspects of phase transitions in s-wave holographic superconductor theory with ordinary 
matter sector coupled to the dark matter one have been elaborated. One expects that dark 
matter authorizes a part of a larger particle sector interacting with the visible matter and 
not completely decoupled [50]-[55]. 

The main point of our studies is the question how the dark matter sector modifies 
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the ordinary phase transitions known from the previous studies of p-wave superconductors. 
The key role will be played by the dark matter coupling constant a, binding dark matter 
helds with the ordinary Maxwell gauge field. It is important to know how the phase 
transitions are modihed by the dark matter sector. To this end the holographic p-wave 
metal/superconductor and holographic insulator/p-wave superconductor phase transitions 
as well as p-wave holographic droplet embedded in magnetic field have been studied. As 
discussed in detail in section 2 the holographic p-wave superconductors provide a numerous 
non-trivial possibilities to look for the influence of the dark matter on them. Out of three 
possible models of holographic p-wave superconductors we shall discuss in detail two of 
them. 

The paper is organized as follows. In section 2 we describe two of the three possible 
models used to build the p-wave symmetry superconductors and conclude that the effect of 
dark matter in one of them is formally the same as for s-wave symmetry. This prompts us to 
consider the other model for which this influence is much more interesting. In the subsection 
B we discuss the crucial points of the considered p-wave holographic superconductor model 
with the influence of the dark matter sector. In section 3, as a gravity background we assume 
hve-dimensional AdS Schwarzschild line element and study the metal/superconductor phase 
transition. Dark matter sector effects on the insulator/holographic p-wave superconductor 
are analyzed in section 4, while section 5 is devoted to p-wave holographic droplet in the 
presence of a constant magnetic field. We conclude our researches in section 6, paying 
attention to the new features of the elaborated phenomena induced by the presence of 
dark matter sector which can potentially serve as an indicator for the future experiments 
dedicated to the detection of dark matter. 

2 Models of p-wave holographic superconductor with dark matter sector 

In the literature on the subject there exist at least three different possible ways of building 
p-wave holographic superconductor. Namely, the Maxwell vector model [23], the SU(2) 
Yang-Mills [11, 15] one and the helical p-wave model [56, 57]. The contemporary review 
of the subject can be found in [58]. Contrary to this, the quantum held theory approach 
(the weak coupling) delivers the unique description of p-wave superconductors, under the 
condition that one ought only to preserve the required symmetry. Namely, from the fact that 
superconductivity is related to a pairing of the two fermions it follows that the total wave 
function of the Cooper pair has to be antisymmetric, with respect to their exchanges. The 
antisymmetry of the spin part requires a symmetric orbital part of the wave function and 
one ends up with s-wave or d-wave superconductors. On the other hand, p-wave symmetry 
of the orbital part of the wave function requires the triplet character of the spin part. On the 
technical level, the held theoretical description of all types of superconductors stems from 
the same type of BCS-like equations which only differ by the symmetry of the form-factors, 
^(k) = 1 for s-wave and ^(k) = kx, where k is a wave-vector, for the simplest p-wave 
symmetry. At present, it is not clear which of the aforementioned holographic models is 
the proper one for the description of strongly coupled p-wave superconductors and what 
are the differences between their properties. 
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In this paper we shall study the Maxwell vector and the SU(2) Yang-Mills models. Our 
analysis relies on the theory in which the gravitational action is given by 

Sg = j d^x [r - 2A^ , (2.1) 

where A = 6/L^ stands for the cosmological constant, while L is the radius of the considered 
AdS spacetime. 

Before we proceed to the main subject of the paper let us give some remarks about 
the other model of p-wave holographic superconductor, the so-called Maxwell vector model 
presented in [27], supplemented by the dark matter sector. The form of the action is in 
a close resemblance of the quantum electrodynamical p-meson, without irrelevant neutral 
part of it [59]. 

The gravitational part is the same as presented earlier, while the matter sector will be 
provided by the action 

Sm = I J (2.2) 

pU pI + i qio Pa pl , 

where a complex vector field with mass m and the charge q was introduced, p^i/ is 
defined by the covariant derivative — iqAfj_ in the form given by 

Pfit^ ~ B)DpPg^. ( 2 - 3 ) 

The last term in equation (2.2) describes the magnetic moment of the vector field p^. The 
vector field constitutes a charged C/(l)-gauge field and on the AdS/CFT side is dual to an 
operator carrying the same charge under the symmetry in question. A vacuum expectation 
value of this operator will be subject to the spontaneous U{\) symmetry breaking. The 
condensate of the dual operator breaks the 17(1) symmetry and moreover because of the 
fact that we have to do with vector fields, the rotational symmetry is broken by choosing 
a special direction. In the light of the above claims, the vector field is treated as an order 
parameter and the model mimics p-wave superconductor. 

Further, let us assume that we shall elaborate real vector field to have the connection 
with the results obtained in our previous studies. By the direct calculations it can be 
checked that, if we assume that the condensate picks out x-direction and 17(l)-gauge fields 
have only f-components 

Pa dx°‘ = px dx, Afj, dx^ = (j){r) dt, By dx^ = p(r) dt, (2.4) 

the underlying Maxwell vector p-wave holographic superconductor with dark matter sector 
and with the real components of the vector field give us the same description of the phase 
transitions as the s-wave model studied earlier [49]. All the equations of motion are of the 
same forms when we exchange 'ip{r) (which acts as an order parameter in s-wave case) for 
the p 3 ;-component of the vector field. Due to this fact we shall not elaborate the vector 
model and concentrate on the SU(2) Yang-Mills one in the following sections. 
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2.1 Model of SU(2) Yang-Mills p-wave holographic superconductor with dark 
matter sector 


In this section we shall describe the basic features of the SU(2) Yang-Mills holographic p- 
wave superconductor model with the dark matter sector. To begin with one makes ansatz of 
an SU{2) Yang-Mills field and two U{1) subgroups of the SU (2), considered as the ordinary 
Maxwell one and the other supposed to describe the dark matter sector, coupled to the 
Maxwell electrodynamics. Next, a gauge boson generated by the other SU{2) generator 
and charged under 17(1) Maxwell subgroup will be taken into account. On the other hand, 
the matter sector is provided by the action 


Sm = I V^d^x(^- , (2.5) 

where and are two SU{2) Yang-Mills field strengths of the form = 

-I- The totally antisymmetric tensor is set as = 1. The 

components of the gauge fields are bounded with the three generators of the SU (2) algebra 
by the relations A = A^^'^ r“ dx^, where [r“, r^] = 6“^^^ The parameter a describes the 
coupling between ordinary and dark matter U (l)-gauge fields. 

The equations of motion imply 


y I “y pallia) , abc p (b) p/iuic) , “ aftc p (b) p/J-uic) _ q 

M 2 M ' 2 ^ 


( 2 . 6 ) 


and for F^i, are provided by 

^abc j^^{b) « pibc ^^{b) p^iu{c) ^ (2 7 ) 


In order to simplify the above equations we multiply relation (2.6) by a/2 and extract the 
term The second term in the equation (2.7) is replaced by the aforementioned 

outcome. The hnal result may be written as 

2 

a V B B (2 8) 

C2 1-^ ^ ^ 

pabc ^ (b) pav{c) “ afec ^ (b) — Q 

I ^fi I 2 ’ 

where d = 1 — 

Both SU{2) Yang-Mills helds, A^^^^ and Bpd>)^ are dual to some current operators in 
the four-dimensional boundary held theory. In order to accomplish a p-wave superconductor 
with dark matter sector we postulate that the following is satished 


A = (p{r) dt + w{r) dx, (2-9) 

B = r]{r) dt. (2.10) 


In the above relations the 17(1) subgroups of SU{2) group generated by are identihed 
with the electromagnetic I7(l)-gauge held {(j){r)) and the other 17(1) group connected with 
the dark matter sector held {r]{r)) coupled to the Maxwell one. The gauge boson held (rc(r)) 
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having the nonzero component along x-direction is charged under A). = (j){r). According 

to the AdS/CFT dictionary, (?i>(r) is dual to the chemical potential on the boundary, whereas 

w{r) is dual to x-component of a charged vector operator. The condensation of w(r) field 

will spontaneously break the U{1) symmetry and is subject to the superconductor phase 

transition. It breaks the rotational symmetry by making the x-direction a special one 

and inclines the phase transition. The transition in question is interpreted as a p-wave 

superconducting phase transition on the boundary. As far as the [/(l)-gauge field bounded 

(3) 

with the dark matter sector is concerned, it has the component ' dual to a current 
operator on the boundary. 

One can remark that the choice described by the relation (2.9) is the only consistent 
choice of the gauge held components allowing the analytic treatment of the problem. The 
direct calculations reveal that the x(l) and t{3) components of the equation (2.6) are given 


as follows: 

« « gi32 ^^(3) ptx{2) ^ ^2.11) 

I = 0, (2.12) 

while the same components of the equation (2.7) imply 

^ gl6c ^^(b) p^^x(c) ^ ^2.13) 

I = 0, (2.14) 

and the main relation (2.8) reduces to the following: 

d ^ 6^32 ^^(3) p.tx{2) ^ gifec ^^(6) ^ ^2.15) 

a = 0. (2.16) 


We take into account the x(l) and t(3)-components of the aforementioned equation. 


3 P - wave holographic metal/superconductor phase transition 


In order to analytically investigate metal/superconductor phase transition we shall im¬ 
plement the Sturm-Liouville method, which for the hrst time has been successfully used in 
holographic phase transition studies [60]. To commence with, one considers the background 
of five-dimensional black hole given by the line element 

ds'^ = -g{r) dt^ + —r {dx^ + dy"^ + dz^), (3.1) 

g{r) L2 

where g{r) = r"^ jl? — r^jr'^L?. Without loss of generality we set F = 1. 

The Hawking temperature for the black hole in question, is equal to Tbh = The 

t (3) and x (1) components of the equation (2.8) are provided by the following relations: 


dr^{r) + 

dlw{r) + 


- dr(l){r) - 
r 



w'^{r) 
a r"^ g 


(t){r) = 0, 


drw{r) -(- 


4>{r) w{r) 
a g^ 




(3.2) 
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The above set of the differential equations should be completed by the adequate boundary 
conditions. We impose that on the black hole event horizon 0(r+) = 0 and the condens¬ 
ing field has a finite norm. The last condition requires that rr’(r+) should also be finite. 
Consequently, the boundary conditions on the event horizon may be cast in the form 

*'(’■) = (i - v) + ■ ■ ■ • P-®) 

“'('■) = '"Zi +Til (i - + ■ ■ ■ ■ (s-*!) 


On the other hand, on the boundary r ^ oo of the considered spacetime one has 

( 2 ) 

(/)(r) -t h - u!(r) (3.5) 

where /r and p are dual to the chemical potential and charge density, respectively, and 
are dual to the source and expectation value of the boundary vector operator. The 
requirement of having a normalizable solution implies that one sets the source equal 
to zero. 

In z = r_|_/r - coordinates the equations of motion have the forms 


4>"{z) - 
w"{z) + 




'Up‘{z) 
a z'^ go 

w'{z) + 


0(z) = 0, 

r\(l){z)[(t){z) - ^giz)] 
a z^ gl 


w{z) = 0, 


(3.6) 

(3.7) 


where ' means derivative with respect to z-coordinate and go = 1 — z^. 

For T ^ Tc the condensate is very small w{z) —)• 0 and w‘^{z) ~ 0. The value of the 
horizon radius for the black hole with temperature Tc is denoted by r+c. The equation (3.6) 
for the (p field reduces to the relation 





(3.8) 


which has the general solution of the form (l>{z) = ci -\-C 2 Z^. Using the boundary conditions 
for it as described by the equations (3.3) and (3.5) we find (j) ~ prZ^^{l — z'^). We need 
solution for the function w{z) near the boundary z —)• 0 of the considered spacetime. We 
shall seek it in the form 

w{z) | 2 ^o=< ^ > z^F{z), (3.9) 

where we have assumed = 0 and denoted < TZ >= 

To proceed further, we shall consider the equation (2.12) in order to find the relation 
between (j){z) and g{z). It can be easily solved and one gets 

viz)+ + C2, (3.10) 


where ci and C 2 are constants. This equation (for a = 0) shows that the asymptotic 
behavior of the field g[z) has the form analogous to (i){z). This is correct as giz) is the 
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counterpart of the field ^(z) in the dark sector. Due to the formal symmetry between both 
fields we assume the similar boundary conditions for r]{z). First, we require ri{l) = 0. It 
leads to the relation 

viz) = z"^), ( 3 . 11 ) 

Comparing this dependence with the behavior of 4>iz) close to the boundary (z —0), we 
introduce the dark matter density pn = ci/4 and approximate the last equation close to 
z = 0, by the following: 

^( 2 ) =-^(1-4) « ^(1-4). (3.12) 

^+c ^+c 

This allows us to rewrite Eq. (3.10) as 

viz) = ^{1-z^)-^(piz). (3.13) 

4c 2 

Introducing the above relation into equation (3.7) and using representation (3.9) of wiz) in 
terms of F{z), we obtain the equation 

ipiz) f'{ z))' - qiz) F{z) + r{z) F{z) = 0, (3-14) 


where we have set 


Piz) = z^ go, 
qiz) = -z^ go f 4 ^ 


r{z) = 


Z^ z go 
zjl- z"^)^ 

90 


and denoted 


with /3 given by the following expression: 


a'^ PD 
4 p 


a r. 


6 ’ 
+C 


/3 — 1 + 


a 


(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 


The equation is solved by means of the Sturm-Liouville method with the function F{z) 
satisfying the conditions F{0) = 1 and F'{0) = 0. The Sturm-Liouville method enables us 
to find the minimum value of A = A^m by minimizing the functional 


^2 ^ fo dz [F'jz)'^ pjz) + qjz) F^jz)] 
fo dz r{z) F^{z) 

The trial function has been chosen in the standard form F{z) = 1 — a z'^. The dependence 
of gir) on r+ has been taken into account by our definition of goiz) = z~‘^ — z^ in the 
functions piz), qiz) and r(z). 








a 



T/Tc 


Figure 1. (color online) Left panel shows the a dependence of the transition temperature 
Tc{a)/Tc{0) normalized to its value without dark matter. The right panel depicts the dependence of 
the condensation value <TZ > normalized to its value < TZ >o calculated for a = 0 on temperature 
T normalized to the transition temperature Tc, for a few values of the ratio x = pd/P and a = 0.5. 
The inset shows the dependence < TZ > / < TZ >o on T for a; = 0 and a = 0, 0.5, 1.0. 


Having in mind that for a given black hole Tc = h can be verified that the 

critical temperature is given by 


Tc{a) 


1 

P3 


vrA®, 


/« _ oIer 

r ^ ^ 

\ a 


1 
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(3.21) 


In the present model of p-wave metal/superconductor phase transition, we find that Tc 

depends on the a-coupling constant of dark matter sector even in the probe limit. This 

dependence reduces to the (a//3)“^/® form for pn = 0, which implies that the bigger a we 

take into account, the bigger value of the critical temperature one obtains. On the contrary, 

without backreaction, we have no such effect in s-wave holographic metal/superconductor 

phase transition in the theory with C/(l)-gauge dark matter sector [48]. For general values 

of dark matter density the behavior is more complicated and is shown in figure 1, for a few 

1 

values of pn- The value of Tc(0) = — is the transition temperature of the considered 

superconductor in the absence of dark matter. 

The previous experience from the study of dark matter sector effects in the backreacting 
s-wave holographic superconductor background [47], makes us to believe that the effects 
should be much stronger for p-wave symmetry superconductors, but this will be a subject of 
the future studies. Evaluating the Sturm-Liouville functional (3.20) we have found Xmin ~ 
4.091 for a ^ 0.6845. For p = 1 one gets Tc(0) Ri 0.199. In the left panel of figure 1, 
the dependence of the transition temperature on the coupling to the dark matter has been 
shown for a number of values Pd/p- We have normalized Tc by dividing it by Tc{0). 


3.1 Condensation value 

In this subsection, we shall pay attention to the influence of the dark matter sector on the 
condensation operator for the phase transition in question and its temperature dependence 
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for T not far from the critical one, Tc- One assumes a small but finite value of the condensate 
< TZ > operator and writes 

=-^{1-z'^)+ <n> ^{z). (3.22) 

nc 

The function ^{z) corrects the function 4>{z) evaluated at Tc close to 2 ~ 0. At temperature 
T the black hole is characterized by the horizon such that T = r^jir. The general 
behavior of the field 0(z) = -^(1 — z‘^) is again deduced from Eq. (3.5). Expanding ^{z) 
around z = 0 

az) = e(0) + C'(0) z + ^e"(0) + ..., (3.23) 

one deduces from (j){z = 1) = 0 that ^(1) = 0. Comparing the coefficients multiplying the 

same powers of z in the two expressions for (j){z) we get 

4 = -^+ < ^ > ?(0) (3.24) 

4 = -^-i<^> ^"(0) (3-25) 

n 2 

c'(0) = 0 (3.26) 


Consistency of the above equations requires ^(0) = —^"(0)/2. To find temperature depen¬ 
dence of < 7^ > we need the value of ‘^^^(0). 

In order to obtain ^^^(0) we use the equation (3.6), introduce the expression (3.22) into 
it and get the relation for ^{z), valid up to the second order in TZ 


" = F\z)il-z^) 

z a rlrl^ go 


(3.27) 


It can be easily verified, by multiplying the above equation by z, that ^^(0) = 0. It also 
provides that ^^^(0) is finite. To find ^^^(0) we rewrite left-hand side of the last equation as 
d{^'{z)/z)/dz and use the boundary conditions ^^(1) = 0. It implies 


C"(0) 


i'{z) , X z{l-z^) F\z) 

- U-^-o— —w / az - 

2 a Jo g 


P<J^E 

2arlrl, ’ 


(3.28) 


where we set 


E = 


f 


F^z) 


irr^- (3.29) 

(1 + Z-^) 

Making use of Eqs. (3.25) and (3.28) we find the dependence of the condensation value on 
temperature close to Tc 


<TZ>= 




T 



(3.30) 


Using the previously found value of a = 0.6845, we obtain E ~ 0.054046 and temperature 
independent prefactor calculated for a = 0 is < 77. >o~ 3.803. 













The right panel of figure 1 illustrates the dependence of the condensation operator 
< TZ > normalized to its a = 0 value < 7?. >o on temperature T, normalized to the 
actual transition temperature for three values of x = pd/p = 0, 3, 6 and a = 0.5. The 
condensation operator is a decreasing function of x, for a given value of T/Tc{a). The 
dependence of the condensation value on temperature for x = 0 and three values of a = 
0, 0.5, 1 is shown in the inset for x = 0. The value of < TZ{T) > decreases with the 
coupling a for constant ratio T/Tc{a) for all values of x. Its behavior as a function of T is 
shown in the inset to the right panel of figure 1 for three values of a. 

The above equation envisages the fact that <TZ > depends on a-coupling constant of 
the dark matter sector directly via factor and indirectly through Tc{a). It also depends 
on the density /9 _d of the dark matter via Tc- The transition temperature increases with the 
coupling |a| for pu = 0. It features a monotonous increase with a for x = 1 and strongly 
decreases with |q;| for the value polp = Q, close to that inferred from astronomical observa¬ 
tions. On the AdS/CFT side, the operator < TZ > can be interpreted as responsible for the 
pairing mechanism. The smaller vacuum expectation value it has, the harder condensation 
happens. So we conclude that dark matter sector destructively influences the condensation 
phenomena in p-wave superconductors for po /p > 1. As visible from figure 1 different 
behaviors are expected for pd/P = 0.1. The opposite (to that in the case po/P > 1) conclu¬ 
sion was achieved in the case of s-wave holographic metal/superconductor phase transition 
[49]. Interestingly, in the Maxwell p-wave model we expect similar dependencies on dark 
matter sector as in s-wave superconductors. The observation of the effect in the laboratory 
could shed some light on the meaning of various models of holographic superconductors 
and their applicability for the description of real life materials. In the case of Gauss-Bonnet 
p-wave metal/superconductor phase transitions [61], the curvature corrections also imply 
the growth of the value of the condensation operator. This fact was also conhrmed by the 
previous studies [62]-[64]. 

4 Insulator/ holographic p-wave superconductor phase transition 

The main ingredient in the considerations of insulator/superconductor phase transition will 
be the gravitational background of hve-dimensional AdS soliton spacetime, given by 

ds^ = dt^ + f{r) dip^ + r^ {dx^ + dy^), (4.1) 

f{r) 

where /(r) = — Xq/x^, xq denotes the tip of the line element which constitutes a conical 

singularity of the considered solution. Without loss of generality we set the radius of the AdS 
spacetime equal to one. The AdS solitonic solution can be gained from the five-dimensional 
Schwarzschild-AdS black hole line element by making two Wick rotations. On the other 
hand, a conical singularity, at the tip x = xq, can be get rid of by the Scherk-Schwarz 
transformation of <y9-coordinate in the form p ~ ip + Ti/rQ. The gravitational background in 
question delivers the description of a three-dimensional field theory with a mass gap which 
in turn takes after an insulator in the condensed matter physics. The temperature of the 
aforementioned background equals to zero. 
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In order to solve the underlying equations of motion for the p-wave holographic in¬ 
sulator/superconductor phase transition problem, we shall impose the adequate boundary 
conditions. Namely, at the tip one has the same boundary conditions as in the s-wave 
insulator/ superconductor problem 

w = wq + wi (r - ro) + W 2 {r - + ..., (4.2) 

(/ = (/o + (/i (r - ro) + 02 (J’ - rof + ..., (4.3) 


where Wa and (j)a, for o = 0, 1, 2,... are constants. One encumbers the Neumann- 
like boundary condition to obtain every physical quantity finite [18]. Contrary, near the 
boundary where r —?> oo, we have the different asymptotical behavior (comparing to the 
s-wave case). The asymptotic solutions read 

w^wo + ^, (4.4) 


where // and p are interpreted as the chemical potential and the charge density in the dual 
theory, respectively. On the other hand, wq and W 2 have interpretations as a source and 
the expectation value of the dual operator. In order to gain the normalizable solution, one 
puts Wo = 0 (we are not interested in the case when the dual operator is sourced). 

In ^-coordinates (with tq = 1) the equations in question yield 


w"iz) + 

r{z) + 


7]{z)] 


+ — 5 /(.) 
f'{z) 1 \ , w‘^{z) 

/(i) +1) ^ = “■ 


w{z) = 0, 


(4.5) 

(4.6) 


4.1 Critical chemical potential 

It was revealed numerically [18] that when the chemical potential exceeds its critical value 
the condensation operator turns on. In the case when p < pc the scalar field w{r) is equal 
to zero, which is interpreted as the insulator phase. The system under consideration has 
a mass gap bounded with the confinement in (2 -|- l)-gauge theory via the Scherk-Schwarz 
compactification. It justifies the conclusion that the critical value of the allowed potential is 
the turning point of the insulator/superconductor phase transition. Moreover, when p ~ pc, 
the scalar field w{r) is so small that w‘^{r) = 0 and the equation (4.6) has an analytical 
solution in the form of a logarithmic function (j)(z) = Ci-|-C' 2 [log(l-)-- 2 ^) — log(l — z^)] [23]. 
Establishing the value of the integration constant at the tip z = 1 (C 2 = 0), one gains that 
4>{z) has the constant value Ci = p. From the boundary conditions, it can also be seen 
that near 2 = 0 we get p = 0. 

As in the preceding section we find the dependence of (p{z) on the component of the 
dark matter sector r]{z). Using the the adequate components of the metric tensor for the 
line element (4.1) and equation (2.12), we arrive at the following relation 


viz) + I (piz) = Di 


z‘^ 1 

y + jlog 


1 + Z^ 

1^ 


+ D 2 , 


(4.7) 
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where Di and D 2 are integration constants. Setting the integration constant Di equal to 
zero at the tip z = 1, we obtain 

D 2 = fJ-D + — (4.8) 

where denotes the chemical potential of the dark matter. For /.i tending to its critical 
values /ic from above, the equation of motion for w{z) function is of the form 


w 


\z) + 




- § Me)] 




a 


w{z) = 0. 


(4.9) 


To proceed we assume wq = and correct the solution for w{z) close to the boundary z = 0 
by dehning the trial function G{z) 


w{z) ~< O > z^ G{z), 


(4.10) 


with the boundary conditions G(0) = 1 and G'{0) = 0 and < O >= W 2 . For simplicity we 
set ro = 1. It implies further, that the equation (4.9) can be rewritten as 

(Piz) G'iz))' - Q{z) G{z) + R{z) F{z) = 0, (4.11) 


where we have dehned the quantities as follows: 


ip 


(4.12) 


a V Me J 


Piz) 

= z^ f, 

(4.13) 

Q{z) 


(4.14) 

Riz) 

= ^3. 

(4.15) 


The minimum eigenvalue of can be achieved from the variation of the functional 


2 ^ fo dz [G^z)"^ P{z) + Q{z) G‘^{z)\ 
fo dz R{z) G2(z) 


(4.16) 


With the choice of G{z) = 1 — a as the trial function we hnd Mmin ~ 2.2666 for a ~ 
0.3376. Let’s note that for a = 0 we have Me ~ should be. The dependence 

of Me on a requires knowledge of which is unknown. For illustrational purposes we 
shall calculate the dependence ^(a) for two values of /i/j = 0, Me- In the hrst case we 

end up with the dependence m — Mmin\J ot/P while in the second case it has the form 
M = Mmin\Ja/[P - a2/4]. 

The dependence of the critical chemical potential on a is shown in the left panel of 
figure 2, for both values ol md- For = 0 the chemical potential is the decreasing function 
of a, independently of its sign. Contrary to that for MD = Me the critical chemical potential 
increases for negative values of the coupling to dark matter, while decreases for positive a. 
Thus, the transition appears at lower values of the chemical potential /ic for positive values 
of a, independently '\i md vanishes or equals that of normal matter. This confirms the fact 
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that dark matter makes the p-wave phase transition easier to happen provided the coupling 
is positive. Let’s recall that for s-wave symmetry the holographic insulator/superconductor 
phase transition the chemical potential was unaffected by the presence of dark matter sector 
[49]. The result was similar as in Einstein-Maxwell theory [18, 24]. 

On the contrary, in the Gauss-Bonnet gravity the critical potential increases with the 
growth of the curvature corrections [26] , making the condensation harder to form. This phe¬ 
nomenon was also found in s-wave holographic insulator/superconductor phase transition 
in Gauss-Bonnet system. 


4.2 Critical phenomena 


When /r —)■ /Xc the condensate value < O > is small but hnite and the equation of motion 
for the (p field is provided by 




/' 1 
— + - 

f z 


< O >2 




a f 


= 0 . 


(4.17) 


For fi slightly above the critical value the condensation scalar operator < O > is very small. 
This enables us to look for the solution in the form 


(j){z) ~ Me + <0 > xiz) + • • • 


(4.18) 


Moreover, to recover the previous result (f){z) = fi, one imposes the boundary condition 
x(l) = 0 at the tip of the considered soliton. Near the boundary z = 0, one expands the 
function x(z) = t(0) + x'i0)z + \x"{^)z^ + ■ ■ ■, and rewrites (j) in the form given by 

(j){z) ~ M - P - Pc+ <0> (^x(O) x'{0)z + ^x"(0) z'^ + (4.19) 


where we have imposed the boundary conditions at the tip of the soliton for the function 
x{z) requiring that x(l) = 0- Comparing the coefficients of the z^-terms, in the above 
equation, we hnd that 

<0>X(0). (4.20) 

Having in mind the relations (4.18) and (4.17), one can easily find the relation for xiz)- 
Namely, it is given by 


X"(z) + ( J + )) X'{z) - 


<0> iicZ^ F‘^{z) 
a f 


+ 0{< O ^ 0 , 


(4.21) 


where by 0{< O >”-2) we denoted terms of order n > 2 which can be neglected because 
of the fact that they are significantly smaller than the linear one. By virtue of the above 
we arrive at the relation 


x"(2) + (J + )) x'iz) - 


<0> Me i^H-2) 


a f 


~ 0 . 


Let us redefine x{z) function by the new one ^(z) 


X{z) = -r- i{z). 


a 


(4.22) 


(4.23) 
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The new definition of x(^) allows to get rid of the d-coupling dependence in the relation 
(4.22). One obtains 




.. F\z) 


Consequently, the scalar operator < O > provides the following: 


(4.24) 


where ^(0) is given by 


<o>= 


^(0) = ai- ^ • 


'(/X - He) a 

f^cm ’ 


(4.25) 


(4.26) 


The integration constants ai and 02 , are determined by the boundary conditions imposed 
on x( 2 )-function. To hnd them we have to fulfill the requirement ^(1) = 0. This leads to 
02 = 0 required to cancel the 02 log(l — 2 ^) term resulting from the integration over 2 ; and 
oi = (3 — 2a + 2a^ — 4alog 2)/24. Using the previously found optimal value of a = 0.3375 
one ends up with ^(0) = ai ~ 0.0894. 

On the other hand, the above relations reveal that the operator < O > yields 


< O >~ A (/X - /Xc)2, 


(4.27) 


where the prefactor A = ~ 3, 139{q:[/ 3 — (a^/4)(/X£)//x)]}^/^ contains informa¬ 

tion on the dependence on dark matter sector. The previously found relation /Xc(d) = 

/x^/d/[/3 — (a^/4)(/X£)//x)] has to be taken into account. 

Equation (4.27) envisages that the p-wave holographic/superconductor phase transition 
represents the second order phase transition for which the critical exponent has the mean 
field value 1/2. This was also the case in s-wave insulator/superconductor phase transition 
influenced by the dark matter sector [48] . Thus the dark matter sector does not change the 
order of the transition. 

The right panel of hgure 2 illustrates the dependence of the condensation value < O > 
on /X close to Hc- The dependence on h/Hc{<^) shows that the amplitude of < O > diminishes 
with a. In fact the dark matter increases the value of the condensate operator < O > if /x 
is kept constant as is visible from the inset. This conclusion is independent of the value of 
Hd provided a > 0. 

Next, we proceed to hnd the dependence of the charge density p on the critical chemical 
potential. In order to calculate charge density p, we consider z^-order coefficients in the 
relation (4.19), having in mind that ■^^(0) = 0, together with the previous requirement 
^(1) = 0 being subject to the boundary condition. It leads to 

y'(0) = ^ 1 ,^ 0 = - < O > ^ fdzz^F\z). (4.28) 

2 ; a Jq 

The comparison of the adequate coefficients of 2 ;^-order in equation (4.19) leads to the 
conclusion that the charge density implies 

P = x''(0). (4,29) 
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Figure 2. (color online) Left panel shows the dependence of the critical value of the chemical 
potential measured in units of jl, for the insulator - holographic p-wave superconductor transition 
vs. a for two values of the dark matter chemical potential = 0, ^c- The dependence of the 
condensation value < O > normalized to its a = 0 amplitude Aq on the normalized chemical 
potential ^cip/) is shown for fijj = 0 and a few values of the coupling to the dark matter sector 
a = 0, 0.5,1.0. The inset to the right panel shows similar dependence on the bar value of and 
for /i£) = 0. 


By virtue of the equations (4.28) and (4.29), having in mind the relation (4.25), one arrives 
at the following 

p={ti- Pc) £>, (4.30) 

where the quantity D yields 


D = 


2 {( 0 ) 


f 


dz F‘^{z). 


(4.31) 


Evaluating the integral with the previously found parameters and introducing ^(0), we 
obtain D ~ 0.844. 


5 P-wave holographic droplet 

In this section, based on the notion of marginally stable modes of vector perturbations 
we shall investigate the influence of the magnetic held on holographic p-wave insula¬ 
tor/superconductor phase transition. We restrict the considerations to the probe limit. 

To begin with, let us recall that the stability of a spacetime can be studied by the 
quasinormal modes (QNMs) of the perturbations in the background in question. When 
one encounters that the imaginary part of QNMs is negative, then they will decrease with 
the passage of time and disappear. The conclusion is that the elaborated spacetime is 
stable against these perturbations. On the other hand, the positivity of the imaginary part 
of QNMs reveals that the spacetime is unstable against the perturbations. The signal of 
instability or possible phase transition in the spacetime in question, is the occurrence of 
marginally stable modes (a; = 0). They are only dependent on radial coordinates and do 
not backreact on the other helds. 
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As the gravitational background we shall use the line element of the five-dimensional 
AdS soliton, rewritten in the coordinates (t, r, p, p, 9). It has the form 


ds"^ = dt^ + 


dr'^ 

f{r) 


+ fir) dip^ {dp^ + dO"^), 


(5.1) 


with /(r) = — r^/r^. We assume the following components of the gauge helds A® and 


B, 


ib). 


A = {pc dt + -Bf p^ dO) P, p) dp, 

B = pd dt, 


(5.2) 

(5.3) 


where Bf is the constant magnetic field. The time components of both helds pc and po 
have constant values. This ansatz is valid close to the critical point of the phase transition. 
As in a previous section pc is the critical chemical potential for the phase transition in a 
droplet, and pD plays the role of the chemical potential of a dark matter. 

The /9(1) component of the underlying equations of motion is provided by 


./, I 1 o2„/, 


1 


dri’+y— + -+ + f 


a 

pc ( Pc ^ pD 


Bj p^ 


1^ = 0. (5.4) 


In order to solve the above equation we choose an ansatz for held which implies 

^ = F{r, t) Hip) U{p). (5.5) 

We set U (p) = p = 2f -^yBf, because of the fact that one looks for the solution conhned to 
a hnite circular region which radius is proportional to 1/ y^Bj. Consequently, we arrive at 
the relations 




drF - - diF + 


1 f Pc jpc - X d-p) Bf 


^2 f 


a 


f 


a 


H 


F = (b.6) 


(5.7) 


One considers the case when A = 2rore/L, where n £ Z causing the periodicity in Hip) = 
Hip + TVLfro). In what follows without loss of the generality we set tq = 1 and L = 1. 
Just, from periodicity property of Hip) we identify A = 2n. 

In the next step, one substitutes T(r, t) = Bir). Consequently, the requirement 
concerning marginally stable modes leads to the condition that cu = 0. Next, the redehnition 
of the r-coordinates in terms of 2 : ones, provides the relation for Riz) 


dzf , 1 


d^Riz) + ( — + -j d^Riz) + ^2 f 


1 


Pc iPc - X ^^d) - Bf Arp 


a 


f 


Riz) = 0. (5.8) 


Then, one introduces a trial function in the form 

Riz) l 2 -^o~ <0 > z'^ 0 ( 2 ), 
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(5.9) 














with the boundary conditions imposed on ©(z). Namely, 0(0) =1 and 0'(O) = 0. 

After some algebra, the resulting equation can be converted into the standard Sturm- 
Liouville eigenvalue equation, which can be written as 


dz 



bn{z) Q + dl c{z) 0 = 0, 


(5.10) 


where we have denoted by (5^ the relation 


e2 _ /^c ifJ'c 4 I^d) Bf 


a 


(5.11) 


The subscript n is used to remind that the function h{z) depends on the parameter n. 
The remaining quantities are defined by 


a{z) = f 

bn{z) = -f z'^ 
c{z) = z^. 


+ 2 


dzf 1 
/ ^ 


4 

z3 f2 


(5.12) 

(5.13) 

(5.14) 


The eigenvalues of 6^ can be found by the method of minimizing the functional in question 


■2 _ {pic - X ~ _ lo {&'{z)'^ a{z) + bniz) e{z)^) 

a dz ciz) 02 (z) 


(5.15) 


To minimize the above functional we have chosen as the trial function B{z) = 1 — az'^. 

The Sturm - Liouville minimization gives a ~ 0.942 and 6'^ ~ 3.8279 for n = 1, while one 

gets a ~ 1.0083 and 6'^ ~ 5.6034 for n = 2. It can be seen that not only magnetic field 

influences the condensation but also a coupling constant of the dark matter sector as well 

(3) 

as the component of the dark matter gauge field , i.e., its chemical potential. 

It can be noted that for the dependence on a cancels out. Having found (5^ we 

get the formula for the critical value of the chemical potential for metal-superconductor 
phase transition at temperature zero. Denoting the ratio of chemical potentials by x = 
AiD/lic, one gets 


Me — 



a 

1 — X jA 


1/2 


(5.16) 


The dependence of the chemical potential /ic normalized by its a = 0 value equal to 
5'^ + Bj on the coupling constant a is shown in the figure 3, for three values of the 
dark matter chemical potential. For illustrational purposes we have chosen //£> = 0, jj,c and 
3fic- For = 0 the critical value of the chemical potential /ic is a decreasing function of 
\a\. It takes on constant value for x = 1 and is strongly increasing function of |a| for x > 1. 


6 Summary and conclusions 

In the paper we have analytically investigated phase transitions towards p-wave holographic 
superconductors in which dar matter sector can imprints its presence. All the calculations 
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Figure 3. (color online). The dependence of the chemical potential pc normalized to + ^/ 

on the coupling constant a for the p-wave droplet for a few values of the parameter x = [>-c- 


were conducted in the probe limit, when the gravitational background was described by 
the five-dimensional AdS Schwarzschild black hole metric and AdS solitonic line element, 
respectively. 

It was revealed that by a special choice of the real vector field, the vector Maxwell model 
of p-wave holographic superconductor is equivalent to s-wave description elaborated in [49]. 
All the properties of the considered phase transitions, like insulator/ superconductor and 
metal/superconductor phase transitions, in both models look identically. Of course, one 
should have in mind that the s-wave order parameter is replaced by component of the 
vector field, being in p-wave holographic superconductor an order parameter. Due to that 
fact we have considered in detail the SU{2) version of p-wave holographic superconductor. 
Our ansatz for S'[/(2)-gauge fields includes two 17(1) groups, the subgroups of SU{2), one 
describing the ordinary Maxwell field and the other related to the dark matter sector. 

Using the AdS/CFT approach and studying the properties of holographic superconduc¬ 
tors in the background with dark matter sector we rely on the following three convictions. 
Firstly, we accept that the gauge-gravity duality teaches us about strongly coupled super¬ 
conductors which are produced and studied in the laboratories. Some examples of this 
can be found in [65]. Secondly, the AdS/CFT correspondence relies on the AdS space- 
time, while our space time is rather the dS one. The power of the gauge/gravity duality 
does not depend whether the Universe in which the superconductors exist, itself forms the 
AdS or de Sitter spacetime. It serves as a kind of calculus, which enables tackling strong 
coupling problems in a d - dimensional field theory using perturbative approach in d -|- 1 
- dimensional gravity theory. Third, the dark matter existence seem to be obvious from 
astrophysical data as discussed in the Introduction. If it interacts with ordinary matter as 
proposed in the relation (2.5) and quantified by the coupling constant a, so it modifies the 
behavior of the ordinary matter as visible from the equations (3.2). These modifications are 
found to influence the superconducting transition temperature and chemical potential and 
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could, in principle, be observed as the annual changes of the properties of superconductors 
following the expected annual changes in the distribution of the dark matter [66]. 

Studies of the p-wave holographic metal/superconductor phase transitions reveal the 
critical temperature dependence on a - the coupling constant of the dark matter sector and 
other parameters of the dark matter: its chemical potential fiD and its density pd - both 
unknown. So for illustrational purposes we have assumed that they may take arbitrary, 
albeit physically sensible values. The results depend on these values. For example it 
turns out that for po = 0 the greater value of a one sets, the bigger value of the critical 
temperature one receives, c.f. the left panel in figure 1. On the other hand for pn/p > 1 the 
increase of a results in decrease of Tc- This is an interesting result, which could potentially 
be of interest in possible experiments aimed at the detection of dark matter. However, one 
should remark that, the qualitatively similar dependence (as here for pz) = 0) was earlier 
found in the case of the backreacting s-wave holographic superconductor in the theory 
under inspection [47]. It has to be checked to which extend the effect will be modihed if 
the backreaction is taken into account in the model under consideration. The condensation 
value < TZ > also exhibits dependence on the dark matter coupling constant. 

As far as the holographic insulator/p-wave superconductor phase transition is con¬ 
cerned, it was spotted that the critical chemical potential also reveals the dependence on 
a. For a: = 0 one observes the decrease of pc with \a\ as shown in figure 2. In the case of 
Gauss-Bonnet gravitational background the situation is quite different. Namely, the critical 
potential increases with the development of the curvature corrections. 

The determination of the condensation operator leads to the conclusion that p-wave 
holographic insulator/ superconductor phase transition is of the second order. The order 
parameter is proportional to (p — with the mean held like exponent. Both pc and 

the prefactor depend on a. Due to this fact, the condensation value < O > may be 
increasing or decreasing function of a depending how one measures the distance from the 
critical chemical potential. If the distance is measured from the actual a dependent pc 
(main picture in the right panel of hgure 2) < O > is decreasing function of a. On the 
other hand, for a given value of p, the increase of the coupling to dark matter appreciably 
enlarges the value of the condensation operator. On the other hand, the charge density is 
independent on a-corrections except via pc- 

We have also elaborated p-wave holographic droplet in a constant magnetic held. In 
our studies we have used marginally stable modes of vector perturbations. It has been 
shown that the critical chemical potential of the superconducting droplet depends on the 
dark matter coupling a, the magnetic held Bf and dark matter chemical potential pn- For 
Pd = 0, Pc diminishes as |a| grows. This is similar to the s-wave holographic droplet case 
with dark matter sector. For po = pc there is no dependence on ce, while the increase is 
observed for pD > Pe¬ 
lt is worth mentioning that due to the fact that cosmic as well as on Earth experiments 
of detecting dark matter are in progress, the untypical behavior of SU{2) p-wave holographic 
superconductors during phase transitions in question may give a hint to unveil its secret. 
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